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Intermod ulation distortion due to nonlinear elements in transistors is 
analyzed using Volterra series representation. It is shown that this technique 
is well suited for the analysis of transistor distortion where the nonlinearities 
are small but frequency dependent. An ac transistor model incorporating 
four nonlinearities is briefly described. The nonlinear nodal equations of 
the model are successively solved by expressing nodal voltages in terms of 
the Volterra series expansion of the input voltage. Based on this analysis, 
a digital computer program has been developed which computes the second 
and the third harmonic distortion for a given set of input frequencies and 
transistor parameters. The results compare favorably with measured values. 
This method also enables the derivation of closed form ac expressions for a 
simplified model; these expressions show the dependence of distortion on 
frequency, load and source impedances, bias currents and voltages, and 
the parameters of the transistor. The technique is also extended to cascaded 
transistors, and simplified expressions for the overall distortion in terms 
of the distortion and gain of individual transistors are derived. Finally, a 
few pertinent practical applications are discussed. 

I. INTRODUCTION' 

Solid-state long-haul analog communication systems are being de- 
signed for higher frequencies to meet the growth in demand. One of 
the more critical and significant problems facing the system designer 
is intermodulation noise arising from transistor nonlinearities. Thus, 
an analysis of transistor distortion at higher frequencies is a practical 
problem; this paper investigates the transistor distortion using the 
Volterra series as an analysis tool. 

Transistor distortion has been investigated in some detail previously. 
Many authors have considered the exponential nonlinear relation be- 
tween emitter current and emitter- to-base voltage which is important 
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at low currents. 1 • 2>3 - 4 - B The effect of frequency on this nonlinear source 
alone has been reported. 5 Three nonlinearities (exponential, avalanche, 
and h FB at dc) have been examined by Riva, Beneteau and Dallavolta. 
For currents up to 20 mA and frequencies up to 100 kHz, Meyer 7, 
has developed a more accurate and complex model obtaining the non- 
linearities from /i-parameters. However, he takes into account the fre- 
quency dependence by assuming that the ^-parameters can be written 
as h' + juh". Moreover, he does not take into account avalanche dis- 
tortion, nor has he extended the model to higher currents (100 mA) and 
frequencies (20 MHz). The model described here considers four non- 
linearities; they are, exponential, avalanche, h FE , and collector capaci- 
tance nonlinearities. These nonlinearities are superimposed on a linear 
ac equivalent circuit. 10,11 Much of the initial development of the model 
with three nonlinearities was done by Thomas. 

The transistor model is analyzed using a Volterra series representa- 
tion; this series is a generalization of the power series. In a now classic 
report, Wiener applied this analysis technique 12 to find the response 
of a nonlinear device to noise. 13 Bose has carried the theory further. 14 
Following a series of lectures by Wiener, 15 the theoretical framework, 
higher-dimensional transforms, and optimization with Gaussian inputs 
were considered by Brilliant, 16 George, 17 and Chesler, 18 respectively. 
Barrett 19 has treated statistical inputs. The synthesis problem has been 
examined by Van Trees, 20 who also applied the method to phase-locked 
loops. 21 The technique has been extended to discrete systems, 22,23,24 and 
a class of time-variant systems. 24,25 More recently the theory of the 
convergence of the series has been treated. 2 ' 1 This work relies more on 
George's work on the higher-dimensional transform theory. 1 

Even though much work has been done in this area, the Volterra 
series has not found a wide application in solving nonlinear system 
problems due to several reasons; if the rate of convergence is not rapid, 
the higher-degree terms, which are cumbersome to handle, cannot be 
neglected; hence, it cannot conveniently represent gross nonlinearities. 
It is not simple to invert the multidimensional transforms to the time 
domain, and it is not a useful technique to determine the stability of a 
nonlinear differential equation. 

The Volterra series method does, however, offer certain distinct 
advantages in analyzing transistor distortion. Since transistor distortion 
is frequency dependent, the power series is inadequate to characterize 
it; the Volterra series does indeed represent frequency dependent sys- 
tems. The nonlinearities in the transistors under consideration are 
extremely small so that the second- and third-degree terms suffice to 
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characterize them. Since the output corresponding to sinusoidal input 
signals is of interest, there is no need to find the inverse of the higher- 
dimensional transforms; the output can be expressed in terms of the 
transform of the kernel. The higher-dimensional transforms of the 
kernel are complex numbers when s< = j«,- , where s, is the complex 
variable in the transform domain; hence, these kernels can be numeri- 
cally evaluated using the computer (see Section IV). Moreover, for a 
slightly simpler model closed form ac expressions can be derived. Since 
the kernels retain phase information, this approach will be useful for 
the AM-to-PM conversion problem at IF frequencies. Finally, in an 
amplifier two or more transistors are cascaded; the nonlinear behavior 
of such cascaded transistors is a significant problem. The Volterra series 
approach can be easily extended to study such cascaded transistors. 

II. AN INTRODUCTION TO VOLTERRA SERIES REPRESENTATION 

A brief exposition of Volterra series with pertinent reference to the 
problem under consideration is presented below. For further details 
the reader is referred to the references cited. 

Consider a simple memoryless nonlinear system described by the 
following power series; let y(t) be the output and x(t) the input; the 
system is represented by 

y{t) = c,.r(0 + c 2 [x(t)] 2 + c 3 [x(t)]\ (1) 

where c, , c 2 , c 3 are constants. For a time-invariant system with memory 
(capacitors and inductors in an electrical network), the linear term 
Jci.r(/)j is replaced by the convolution integral (x(t) = 0; t < 0) 

Vi(t) = f d(t - t)x(t) dr. (2) 

Iii the transform domain, (2) may be written 

r,(«) - WW. (3) 

This transform domain representation of the system [Cj(s)] has been 
an invaluable aid to the communication engineers since it brings into 
focus the frequency behavior of the system. 

A generalization of the second-degree term, c 2 [x(t)] 2 , is the double 
convolution integral 

!h(0 = f f *(. - t, , t - r 8 ) II x( Ti ) du . (4) 

J U J 1 = 1 
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The output depends on the past values of the input; the above expression 
involves a product of the input with itself, thus representing a quadratic 
system. c 2 (t — r if t— t 2 ) is known as the second-degree Volterra kernel. 
A two-dimensional Laplace transform can be denned for (4) after 
introducing dummy variables U and U . As shown in Appendix A, (4) 
becomes 

2 

Y 2 (s l ,s 2 ) = C 2 (s 1 ,s 2 ) II *(*•■)• ( 5 ) 

i = l 

When two sinusoidal signals at frequencies /„ and f b are applied 
(Appendix A), the output at the harmonic frequency /„ ± f b is given 
by [| C 2 (/ a ± h) | cos (27r(/ a ± f h )t + a *»)]. Since in general C 2 (/„ , / 6 ) 
will not be equal to C 2 (/„ , — /&), different values of distortion at different 
harmonic frequencies are directly reflected in the kernel. Moreover, as 
in the power series case, the 2/ product is less by a factor of two. 

Likewise, the third-degree term [c 3 (z(r)) 3 ] can be generalized to a 
triple convolution integral; 

2/3(0 - f f f C ^ 1 - T » ' * - T " ' ~ r ») II *(r0 dr, . (6) 

Jo •'0 «^0 *■! 

In the transform domain (6) may be written 

3 
Y a ( Sl , s 2 , s 3 ) = CM , s 2 , s 3 ) II *(*<)• (7) 

The magnitude of the signal at the harmonic frequency /„ + f b — /„ 
due to the three fundamental signals at /„ , f b and f e is given by | C 3 (/ a , 
fb , —fc) I- The constants like 1/4 for a '3/ a ' product are the same as 
obtained from the power series approach. 

Later in the paper (in Section IV) the cascade relations in the trans- 
form domain are frequently used; their physical significance is discussed 
in detail in Section VI. (See also Fig. 1.) The cascade formulae and the 
procedure for deriving them are given in Appendix A. 

The second and third harmonic distortion are denned as the second 
and third harmonic power in dBm, respectively, when the fundamental 
power at the output of the transistor is at zero dBm (one milliwatt). 
In the analysis of the model in Section IV, the output voltage is ex- 
pressed in terms of a Volterra series of the input voltage. Thus, the 
kernels Ci(sO, C 2 (s x , s 2 ), and Cafa , s 2 , s 3 ) are the voltage transfer 
ratios; for a given load R L , the second and the third harmonic distortion 
in dBm are given by the following expressions: 
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(8)* 
(9)* 




Z=F(X) 
Fig. 1 — Two cascaded systems. 



III. THE JUNCTION TRANSISTOR NONLINEAR MODEL 

A model is a simple but realistic representation of a physical phe- 
nomenon in terms of measurable parameters such that the phenomenon 
can be analyzed, and controlled if possible. The linear equivalent circuit 
of a transistor is one such example. In reality, several elements of the 
transistor equivalent circuit are not linear but are linearized versions 
of nonlinear functions; they are the first-degree terms of the Taylor's 
series expansion of the nonlinear functions. Hence, a logical way to 
develop the nonlinear model is to consider the second- and third-degree 
terms of the Taylor's series expansion; thus, the emitter resistance 
(exponential nonlinearity), current gain (h FE and avalanche nonlin- 
earity), and the collector capacitance (collector capacitance nonlin- 
earity) have been represented by nonlinear voltage dependent current 
generators whose parameters are higher-degree Taylor's series terms. 
This approach has another advantage in that it is difficult to measure 
the nonlinearities since they are small; but, it is not too difficult to 
measure the overall functions and to curve fit with the known theoretical 

* The factors \ and \ normalize the distortion to 2/ and 3/ products. 
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Fig. 2 — Common-emitter nonlinear equivalent circuit. 

or empirical relations. These nonlinearities are superimposed on the 
linear equivalent circuit (Fig. 2). The nonlinearities are described next. 

3.1 Exponential Nonlinearity 

The emitter current, I E , is related to the emitter voltage, V 2 , by 
the exponential relation 



I E = a\ exp [q j~) - 1 + B, 



(10) 



where K = Boltzmann's constant, 

q = electron charge in coulombs, 

T = Temperature in degrees Kelvin, 

and A and B are constants which depend on the transistor parameters 
(Ref. 27; p. 181, p. 249). An experimental curve of the emitter current 
I E and the emitter-to-base voltage V b . is shown in Fig. 3. This non- 
linearity is expressed as a voltage-dependent current generator by a 
Taylor's series expansion of (10) as follows: 

i. = K(v 2 ) = K> 9 + K*vl + K 3 vl , (11) 

where the Taylor's series coefficients are derived in terms of known 
parameters, the emitter resistance r,, and the emitter bias current I E ; 
i.e., 
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K __ _L . K I 

1 ~ r ei ' Aa ~ 2(r t JI E 



4-: *.-i X 



1 



6 (/*)- (r.,)* 



(12) 



3.2 Avalanche and h FE Nonlinearity 

The collector current is a nonlinear function of the emitter current 
at higher values of current (h FE nonlinearity) and of the collector-to-base 
voltage at higher values of voltage (avalanche nonlinearity)." 7 h FE , 
the ratio of I c to I B , is plotted as a function of collector current I c 
in Fig. 4. It is seen that the following empirical relation 6 matches the 
experimental result (Fig. 4) : 



hp E — 



hp E , 



(13) 



1 + a log 2 



where h FE max is the maximum value of h FE , I c mux is the value of I c 
at which h FB max occurs, and a is a constant. 

The avalanche nonlinearity is due to avalanche multiplication which 
occurs at higher collector-to-base voltage. It is determined from the 
collector characteristic which is a plot of collector current (/ c ) and 
collector-to-emitter voltage (V C e), (Fig- 5). The empirical Miller's 
avalanche multiplication factor is given by 



1 - 



V< 



Vi 



(14) 
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Fig. 3 — Exponential nonlinearity — measured curve. 



998 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1967 




5 6 8 10 



20 30 40 50 60 80 100 

I c IN mA 



200 300 



Fig. 4 — h F E nonlinearity — calculated and measured curves. 



where V C eo is the sustained voltage, and the exponent n is determined 
by experiment. From expressions (13) and (14), the ratio I c VsIe is 
given by 



Ie 



h* 



1 



1 + h FE max + a log 2 



1 - 



V< 



V CBO 



where V CB o = V C Eo/n\/l — a and V CB = Vce • 



(15) 
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Fig. 5 — Avalanche nonlinearity — calculated and measured curves. 
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The ac i c can be expressed in terms of i, and v cb [v 3 — v x ] by a Taylor's 
series expansion of (15). Since i, is a function of emitter voltage v 2 , i e 
is represented by a current generator g(v 2 , v 3 — v x ); for convenience in 
notation it is separated into a linear term g x (v 2 , v a — v x ), a second- 
degree term g 2 (v 2 , v 3 — v x ) and a third-degree term g 3 (v 2 , v 3 — v x ). 
The linear term equals M (<x x K x )v 2 + M x (v 3 — v x ). The second-degree 
term is given by a 2 M K 2 x (v 2 ) 2 + m 2 (v 3 — v x ) 2 + (a x M x )K x v 2 (v 3 — v x ). 
The coefficients a x , a 2 , M x , m 2 , etc., and the third-degree term are 
given in Appendix B. 

3.3 Collector Capacitance Nonlinearity 

The collector capacitance is a nonlinear function of collector-to-base 
voltage (V CB ) since the depletion layer width is a function of V C b • 
The exact functional relationship is determined by plotting the common- 
base imaginary part of h 22 as a function of collector-to-base voltage 
(V C b) as shown in Fig. 6.° It is evident from the figure that C c follows 
the 1/3 voltage law (Ref. 19; Equation 5-96); 

C e = k(V CB y 1/3 . (16) 
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Fig. 6 — Collector capacitance nonlinearity — calculated and measured curves. 
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This nonlinearity is represented as a frequency (differentiation) and 
voltage-dependent current generator as follows: 

it. = y(v 3 — »i) 

= 7l Jt ^ 3 ~ ^ ~*~ y 'di ^ 8 " Vl ^ ^~ y *di ^ 3 ~ Vl ^' ^ 

where ji = C e , and where y 2 and 73 are known from (16). 

The above nonlinear current generators are incorporated in the linear 
equivalent circuit as shown in Fig. 2. The linear equivalent circuit 
parameters are obtained from the equivalent circuit characterization. 
They can, for example, be computed from the /i-parameters at different 
frequencies. In general, the distortion is not a critical function of the 
linear parameters. (Figs. 14 to 17). 

All the nonlinear coefficients (K 2 , a 2 , wi 2 , etc.) are easily obtained 
from a simple computer program. The parameters to be specified along 
with typical values for transistor type A-2436 are listed in Appendix C. 

IV. THE VOLTERRA KERNELS FOR THE NONLINEAR MODEL 

The Volterra series method is applied to the model to compute the 
second and the third harmonic distortion. The voltage at each node is 
a nonlinear frequency-dependent function of the input voltage. Each 
nodal voltage is expressed by a Volterra series expansion of the generator 
voltage; since the nonlinearities are small only three terms are con- 
sidered. The kernels at each node are determined from Kirchoff's current 
equations. 

4.1 Nodal Equations 

The Kirchoff's current law is applied at each node; the currents are 
next expressed in terms of the generator voltage v g , the three nodal 
voltages »i , v 2 , and v 3 , and the known linear and nonlinear parameters. 
The impedances are represented by their transforms and denotes 
that it operates on the voltage across it. The nodal equations are given 
below. 

' o(v„- y,) + (sC 3 ) o 0-3 - p,) = (sCO v v + - (», - v 2 ), (18) 



- (y, — v 2 ) = K{v 2 ) + (sC 2 ) ov 2 + ( — ) (v 2 — v 3 ) - 7(1*3 — v 2 ) 



(19) 
g(v 2 , v 3 - v,), 
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-y(v 3 — v 2 ) + - (v 2 - v 3 ) — g(v 2 , v a - v t ) = (sC 3 ) o (v s — v,) 

where K(v 2 ), y(v 2 — v 3 ) and g(v 2 , v* — v,) are the nonlinear current 
generators. 

4.2 Solution Using Volterra Series 

Since each nodal voltage is to be expressed in terms of three Volterra 
kernels, there are nine unknown Volterra kernels to be determined from 
the three equations. The problem of solving for nine unknowns from 
three equations is resolved by noting that the polynomials x, x 2 and 
# 3 are linearly independent; hence, each degree term is separately 
and successively solved. The linear kernels are first determined; then 
the second-degree kernels are determined in terms of the linear kernels; 
lastly, the third-degree kernels are evaluated in terms of the first- and 
second-degree kernels. 

Let Ai(s), Bi(s), Ci(s) denote the transforms of the linear kernels 
at nodes one, two and three, respectively. From the nodal equations 
(18) to (20), the following vector matrix equation is derived. 



[ 1 




At(«) 





■ = P»W = 


*.(«) 


o , 




koo 



(21) 



where 



PsOO 



- i i 

+s(C«+Ci) + — 

Z,,(s) n 



h"»i 



— sCi — mi 



-aC, 



1 



hsCsH hKi(l -«) +*7i mi— 8Ti 



1 

\-mKi—Syi 



+«6'3 + 



%(■) 



- +701 + 871 



(22) 

Equation (21) is solved by inverting matrix P E (s) and post-multiply- 
ing by the vector 

J_ 
ZM 



I 
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For a given frequency s = ju>, the computation is done numerically. 
The second-degree terms are equated next in (18) to (20). There are 
two types of second-degree terms; those arising from the unknown 
second-degree kernels [for example, (s, + s 2 )Ci^2(si , s 2 )] and those 
arising from the known nonlinear coefficients and the known linear 
kernels [for example, K 2 n* =1 B,(s,)]. The terms associated with the un- 
known second-degree kernels are the same as were associated with the 
unknown linear kernels in (21), but at the harmonic frequency (sx + s 2 ). 
The following vector matrix equation is obtained for the second-degree 
kernels: 



+Y2(C\ - B t ) - K 2 {B X )] 
-y 2 (Ci - £i)]J 



r = Pb(si + s 2 ) < 



A2{Si , s 2 ) 
B 2 (s l , s 2 ) 

C 2 (Si , s 2 ) 



(23) 



where § 2 and f 2 represent the second harmonic contribution due to 

g 2 (v 2 , v 3 — Vi) and y 2 (v 3 — v 2 ); hence, 

§ % {B l , C x - A 1 ) = [aA + ajQ t[BM 

i-l 

+ ^^ [BMWM - AM] 

+ B^lCM - AM]] 

+ m 2 II [CM - AM] (24) 

R 2 {B X ) =K 2 flBM (25) 

UCi - BO = K.(«i + s 2 ) 72 II [Ci(«.-) - BM]- (26) 

i = l 

J°£(si + s 2 ) is the matrix P E (s) with s replaced by (sj + s 2 ). 

The vector on the left side of (23) is known. Thus, the unknown kernels 
are determined by inverting the matrix P B (Si + s 2 ) and post-multiplying 
by the vector on the left-hand side of (23). When s t = j<a b , the inversion 
of the matrix and the post multiplication by the vector can be done 
numerically. 
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The procedure for obtaining the third-degree kernels is almost the 
same; the significant difference is that the vector on the left side not 
only contains terms arising from the third-degree nonlinear parameters 
but also includes second-degree coefficients which give rise to third- 
degree terms by the interaction of the first- and the second-degree 
kernels. These interaction terms are denoted by K 23 , § 23 , 723 , respec- 
tively. For example, K 9 (Bi) = K 3 n< =1 2?i(s,-), whereas K 23 = 2/c 2 E 1 (s 1 ) 
B 2 (s 2 , s 3 ) which shows the interaction of the first- and the second- 
degree kernels. The third-degree kernels are derived from the following 
equations: 



[MBy , C x - A t ) + y 3 (C 1 - B,) 

- 2t 8 (B0 + <? 23 + 723 + #2 3 ] 

I- MB* ,d- A0 -7 3 (d -B,) 

— <?23 — 72 3 ]J 



= P*(Si + s 2 + S3) 



-^-3\ S 1 J S 2 > S 3) 

B 3 {s x , s 2 , s 3 ) 

.^3(^1 ) S 2 > S3), 



(27) 



where 3 , <? 2 s are given in Appendix B. 

A computer program has been developed which calculates the kernels 
and the second and the third harmonic distortion. It uses existing pro- 
grams to invert the matrix P E (s). The nonlinear coefficients are com- 
puted from the known and measured parameters. Computed and meas- 
ured results at different currents are given in Fig. 7. The program has 
been extended to common-base and common-collector configurations. 

V. SIMPLIFIED DISTORTION EXPRESSIONS, THEIR PHYSICAL SIGNIFICANCE 
AND COMPARISON WITH EXPERIMENTAL RESULTS 

Another advantage of the Volterra series method is that it permits 
derivation of closed-form expressions for second and third harmonic 
distortion. These equations show the interaction between the various 
nonlinear parameters and the effect of frequency. 

The model includes the base resistance (r 6 ), the emitter resistance 
(r„), the diffusion capacitance (C 2 ), the load (R L ) and the source im- 
pedances Zg(s), and three nonlinearities, namely, exponential, ava- 
lanche, and h FB nonlinearities. In the computer program C be , C h . , 
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Fig. 7 — Comparison of experimental and computed results. 

r e , C e , mi and collector capacitance nonlinearity have been taken into 
account. The expressions given below are for the common-emitter 
configuration. 

5.1 The Second Harmonic Distortion Term 

The second harmonic distortion in dBm (8) is given by 

(r b + Z g (g))-(K, + sC 2 ) + 1 



ilf, 



20 log - 






,(r b + Z.(s))-[KiO- ~ «0 + fiCJ + 1- 



fo-4 






+ m 2 n ^ + 






+ 



fo 



(ZM + rJ-gC. + l 



(28) 



«,# \(Z g (s) + r fc )-(X:i + sC 2 ) + l/_ 
where s, = j<a a , s 2 — ±ju t and s = ju a d= ;co 6 . 

5.2 TVie TVwrd Harmonic Distortion Term 

In the third harmonic distortion term given below, the interaction 
terms due to the first- and the second-degree kernels have not been in- 
cluded mainly to reduce the complexity; in certain cases, they may be 
significant. 

(n + Z g (s))-(sC2 + K i ) + l 



M, 



»*i¥ 



An + ZM<KxQ - «i) + aC) + i_ 
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K 



2 R L sC 2 r b 1 SjSjClrl 

3 ai K, + 3 (a.K,) 2 



% + & x [ri + i^p^ + |Hfc^j _ ,„ 3 n (** + 



s,C 2 rb 



otxlin + ZM(sc») + i] + 2a 2 g 2 



(a 1 K 1 ) 3 (r 6 + Z,(s))-(sC 2 +A' 1 ) + 1 (aO^i 



Oj 



M, 



(«:) J 



OiSx + I 



C. 



3 (K,r_ _ 



(29) 



where s, = jw a , s 2 = ja b , s 3 = ±j<a e and s = s, + s 2 + s 3 and s,s, = 
s,s 2 + s 2 s 3 + S 3 Sj . 

5.3 Physical Interpretation of the Distortion Terms 

The interaction of different nonlinearities and their dependence on 
load impedance, source impedance, bias currents, bias voltage and 
frequency is indeed somewhat complex. However, the closed form ex- 
pressions derived above give a general qualitative picture which will 
be discussed now. 

5.3.1 Effect of Frequency 

It is important to know the effect of frequency on distortion. The 
distortion depends not only on the frequencies of the fundamental tones 
but also on the harmonic frequency of interest. As shown in Fig. 8, 
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M 2 due to a + b product is better than M 2 of a — b product by 10 dB 
with the two tones at 8.32 and 7.266 MHz. These measurements were 
made with the transistor biased at 100 mA, 10V and with R L = 50ft 
and R„ = 50fi. In curve (a) of Fig. 8, the fundamental tone was in- 
creased from 2 MHz to 10.5 MHz and signals at 2/ and 3/ were measured. 
It is seen that both M 2 and M 3 improved with increase in frequency. 
A theoretical explanation on the basis of dominant terms (in this range 
of parameter values) in (28) and (29) is given below. In (28) as well as 
in (29) the terms in brackets are multiplied by a frequency-dependent 
term 

(r 6 + «„)(#! + iC 2 ) + 1 



(r„ + fl.)(tfi(l - «i) + §C 2 ) + 1 

In this range of frequency (s = harmonic frequency), if K t (l — «i) ^ 
| sC 2 | and if | (r b + R„)sC 2 \ > 1 but K x > \ (sC 2 ) |, then the above 
term reduces to K t /sC 2 which decreases with increase in frequency. 
However, the avalanche terms (M 2 , M 3 , etc.) involve the terms sC 2 
[in (28) and (29)] and SiC 2 in the numerator. Thus, if the avalanche 
terms are dominant, as at higher voltages, there should be no net con- 
tribution due to avalanche terms alone. The exponential terms [K 2 /(Ki) 2 
and K z /(K t ) 3 ] are multiplied by the factor 

(r t +fl g )-§(7 2 + 1 



(r*+22,)-(£i(l - «0 + §C*) + 1 

This term is independent of frequency if (sC 2 (r b + R g ) + 1) > 1. 
Thus the above discussion shows that distortion will improve with 
increase in frequency at lower voltages and if I sC 2 (r b + R„) + 1 | < 1. 
To verify this statement, the voltage was increased to 20 volts and the 
input resistance changed to 22512. The plots of M 2 and M 3 with fre- 
quency, as measured, are given in curves labeled (b) in Fig. 8. It is 
seen that M 2 and M 3 do not improve with increase in frequency. The 
small improvement can be attributed to the h FE terms. 

In general, increase in frequency increases distortion; this is especially 
true for the common base configuration. But as shown above, for certain 
ranges of frequency and certain values of source impedance, distortion 
can improve with frequency. 

5.3.2 Effect of Load Resistance, R L 

The load resistance is an external parameter which the circuit designer 
can vary; hence, it is useful to know its effect on distortion. The second 
and the third harmonic terms are multiplied by l/\/R L and 1/R L 
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terms, respectively; it shows that the distortion can be reduced by in- 
creasing R L . However, the avalanche terms SL X , M 2 , and Jft 3 are 
multiplied by the R L term, so that increasing R L will increase the con- 
tribution from the avalanche terms. Thus, an increase in R L may in- 
crease distortion or reduce it due to cancellation. (The contribution from 
the collector capacitance terms also increases with increase in load (Rl).) 
Because of the above interaction, for a given set of parameters and 
input frequencies and the harmonic frequency of interest there exists 
an optimum load R L ; this, of course, can be determined using the com- 
puter program. In Fig. 9, the measured values of M 2 and M 3 at different 
values of R L are plotted; in both cases increasing R L reduces distortion 
until the optimum value is reached and then distortion increases with 
increase in R L . 

5.3.3 Effect of Source Impedance, Z„(s) 

Source impedance is another important external parameter. The 
source impedance affects the exponential nonlinearities K 2 /K\ in (28) 
and IU/Kl in (29) by the factor 

(r h + Z a (s))sC 2 + 1 



(r» + Z g (s))-[sC 2 + Kt(l - a,)] + 1 

At low frequencies, this nonlinearity is reduced by the factor 1/[(1 — a) 
(72 a _|_ r(j )2C, + l]. Thus, an increase in R„ will reduce distortion from 
this source. However, the contribution from other nonlinearities are 
increased by 

1 + KAR, + r b ) 











1 + K x {R t + r 6 )(l - «,: 
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It is seen from this expression that if KxiR,, -f- r b )(l — a r ) is greater 
than 1, the other nonlinearities are not affected by the increase in R g . 
Thus, if the exponential term is dominant, increasing R„ reduces dis- 
tortion at low frequencies. At higher harmonic frequencies if 
| sC 2 (Z„(s) + r b ) [ > 1, the distortion terms are independent of the 
source impedance since the [r b + Z„(s)] term in the numerator and 
denominator cancel. This is well illustrated in the measured results of 
Fig. 10. The second harmonic frequency being 0.7 MHz, | sC 2 (R B + r b ) \ 
is not much greater than one up to R„ = 100ft; hence, the second har- 
monic distortion improves with increase in source resistance up to 
140ft. Further increase in R g does not cause much change in distortion. 
The third harmonic frequency is 17.3 MHz; hence, a change in R a does 
not affect M 3 appreciably. ([ (sC 2 )-(R„ + r b ) \ > 1) 

5.3.4 Effect of Bias Current 

Increase in bias current usually reduces distortion due to the following 
reasons. The increase in emitter bias current reduces the exponential 
terms 



K, 



h ] :imi 



k 3 r i 



(K x y V I J (Ki 

Fig. 11 shows the effect of bias current on h FE terms; a 2 decreases 
with increase in I c by 

1 i / he 
T log U 

J- C \-» C urn; 
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it becomes zero at I c = Ic max/e, and then becomes negative, and 
increases with further increase in I c . The coefficient a 3 decreases with 
bias current I E . Thus, in general, an increase in bias current has the 
effect of reducing both the second and third harmonic distortion (Fig. 7) 
(at least until a 2 = 0). 

5.3.5 Effect of Bias Voltage 

Whereas exponential and h FE terms are functions of bias current, the 
avalanche and collector capacitance nonlinearities are affected by the 
bias voltage. The coefficient Jl? 2 increases with the voltage; but i#i and 
M 3 increase much more rapidly (Fig. 12). (Both the collector capacitance 
nonlinear coefficients y 2 , 73 decrease with the increase in bias voltage.) 
The effects of change in bias voltage are especially pronounced at higher 
load resistance since avalanche (and collector capacitance) terms be- 
come dominant. The third harmonic distortion decreases more with 
the increase in voltage (Fig. 13) than the second harmonic distortion 
does. 
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The physical significance of the closed form expressions has been 
qualitatively discussed. Precise quantitative estimates can and have 
been obtained using the computer program. For example, the effect 
of varying linear parameters by fifty percent of their original values was 
studied. The results show that the distortion does not critically depend 
on the linear parameters (Figs. 14 to 17). The other transistor parameters 
such as I c mB x , Vcbo , n, etc., can also be varied. 

VI. ANALYSIS OF CASCADED TRANSISTORS 

It is often stated that in a multi-stage amplifier, the output stage 
alone determines the over-all distortion. Even though this statement 
is true to some extent, it is frequently found in practice that the effects 
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of the previous stages cannot be ignored and sometimes the previous 
stage is dominant. This is especially true if both minimum noise figure 
(which requires lower bias current) and modulation requirements are 
to be met by a two stage amplifier. Two analysis tools based on Volterra 
series are presented here which enable the study of such cascaded stages. 

The first approach makes use of the cascaded formulae mentioned 
earlier; this method illustrates the cascade phenomenon and permits 
derivation of simple cascade rules. 

Consider two cascaded transistors (Fig. 1); let the output voltage 
(v 2 ) of the first transistor be denoted by D(v,); the output voltage (v 3 ) 
of the second stage by E(v 2 ) and F(v g ). The aim is to compute the kernels 
F i(«i), F 2 (s 1 , s 2 ), F 3 (s l , s 2 , s 3 ) knowing D and E. To calculate D(s,), 
etc., it is necessary to know the load impedance of the first stage which 
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is the input impedance of the second transistor. This can be computed; 
thus, for a given generator impedance and bias conditions, D^), 
D 2 (s! , s 2 ), D^ , s 2 , s 3 ) can be determined. Efa), E 2 (si , s 2 ) and # 3 (si , 
s 2 , s 3 ) can be computed for a given load and bias conditions with R„ = 
(voltage v 2 is directly impressed across the second). Now expression 
v 3 in terms of v„ is given by 



v 3 = F(v.) = E(v 2 ) = E(D(v,)) = {Eo D)(v,). 



(30) 



It is seen that F is related to E and D by the cascade formulae whose 
physical significance is discussed below. 



6.1 Linear Term 

The linear term is given by 

FM = D t (s)E 1 (s) 



(31) 



which states that the overall gain in dB is the gain of the first stage 
in dB plus the gain of the last stage in dB. 

6.2 Second Harmonic Term 

The second-degree kernel is given by 

2 
F 2 ( Sl , s 2 ) = E^s, + s 2 )D 2 ( Sl , s 2 ) + E 2 (s t , s 2 ) II DM- (32) 



The first term of the formula states that a given harmonic product 
from the first transistor D 2 (jo) a ± ju b ) is amplified by the second transis- 
tor at the harmonic frequency E x (ioi a ± ju b ). The second term shows 
that the two fundamental tones are amplified by the first transistor 
[Di(ji0 a )Di(zkju b )] and then the second transistor acts on these tones 
to produce distortion E 2 (j<j) a , ±;w 6 ). 

Equation (32) is related to the second harmonic distortion (M 2 ) 
as follows: 



M 2 = 20 log 
= 20 log 



V10-X 

2 


F 2 (si , s 2 ) 
FMFM 


VlO~ 3 R L 




-D 2 (s, , s 2 ) 


2 


n DM 

1=1 



EM + s 2 ) EMbi , s 2 ) 

VMEM + JT^ 



(33) 

(34) 



The second term is the second harmonic distortion of the last stage. 
The first term expresses the contribution from the first stage; it approxi- 
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mately equals 

I" First stage second harmonic! — ["gain of the last stage! (35) 
[distortion in dBm J |_ m dB J 

This shows that if the gain of the last stage is high, the contribution 
from the first stage is small. Equation (35) is approximate in two re- 
spects; it neglects the frequency effects and the phase addition of the 
contributions from the first and the second stage. In (35), the second 
stage gain in question is actually the ratio 

#i(Si + a a ) 

E l (a l )EM 

which involves the two fundamental and the harmonic frequencies. 
As an example, a shaping network which was introduced increased the 
gain (18 dB) at the harmonic frequency (0.7) MHz) and decreased the 
gain at fundamental tones 15.2 MHz (8 dB) and 14.5 MHz (8 dB) 
with the result the overall distortion was poorer by 34 dB. 

6.3 Third Harmonic Distortion Term 

The third harmonic kernel F 3 (si , s 2 , s 3 ) is given by 

F 3 (s, , s 2 , s„) = S,(s, + s 2 + s 3 )£ 3 (s, , s 2 , s 3 ) + 2J0 3 (s, , s 2 + s 3 ) 

3 

•Di(8x)D 2 (s 2 , s 3 ) + E 3 ( Sl , s 2 , s 3 ) II DM> (36) 

The first term shows that the third harmonic product of the first stage 
tf>»(8i , s 2 , s 3 )] is amplified by the last stage at the harmonic frequency 
[#i(si + s 2 + 8 a )]. The second term is the interaction term; it arises 
when the second-degree kernel of the last stage [E 2 (si , s a + s 3 )] acts 
on the sum of the fundamental [D^sJ] and the second harmonic output 
of the first stage [D 2 (s 2 , s 3 )]. The last term shows that the second stage 
third-degree kernel [E 3 (s t , s 2 , s 3 )] acts on the fundamental tones am- 
plified at the respective frequencies by the first stage [Di(s 1 )D 1 (s 2 )Di(s 3 )]. 
From (36), the overall third harmonic distortion is related to that 
of the individual transistors by 



M 3 = 20 log | IQ~ 3 R L 



Jffife + s 2 + s 3 ) 



+ 2 



gajgi | S 2 + 8 3 ) 



E 1 (s 3 ) 



Hem 

Q 2 («2 , S 3 ) 

II DM 



D3\ S 1 1 S 2 ) S 3) 



+ 



'II DM 

i=\ J 

E 3 {s x , s 2 , s 3 ) 
fl EM 



(37) 
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The first term is the contribution from the third harmonic term of 
the first stage; it is given approximately by 



Third harmonic distortion 
|_of the first stage in dBm J 



— 2 ["Gain of the last 
L stage in dB 



(38) 



The interaction term approximately equals 



Second harmonic distortion 
[_of the first stage in dBm J 



+ ["Second harmonic distortion ~| 
_of the second stage in dBm J 



+ 6dB - 



Gain of the last 
|_ stage in dB 



1 



(39) 



The third term in (37) is the third harmonic distortion of the last stage 
in dBm. 

It is seen that the effect of the first stage and the interaction term 
can be reduced by increasing the gain of the last stage. Equation (39) 
illustrates that the second harmonic distortion of each stage should be 
good. This may become a limitation if the first stage is biased at lower 
currents. 

In the above simplified expressions [(38) and (39)] phase addition and 
frequency effects have not been considered. In (38), 2 (gain in dB) 
actually represents 



20 log 



■Ei(jco ± ju b ± p) e ) 



In (39) the second harmonic distortion is to be measured with two 
tones, one at the fundamental and the other at the harmonic frequency 

-E 2 (Si , S 2 + 8 a ) 
E l (8 l )E l (8 2 + «,) 

and then multiplied by the ratio of the gain 

#l(*2 + s 3 ) 

Moreover, the kernel must be made symmetrical by taking the average of 
three possible combinations. 

Thus, the simplified expressions (35), (38), and (39) are exact if the 
transistor performance is not frequency dependent; in general, they 



1016 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 19G7 

can be used to get a qualitative picture. Equations (34) and (37) are 
indeed exact and take into account the frequency dependence. The 
computer program is being extended to calculate (34) and (37). 

An alternate approach to calculate the distortion of cascaded stages 
is to analyze the nonlinear equivalent circuit of cascaded transistors 
using the nodal technique illustrated in Section IV. The nodal equations 
are derived first; next each nodal voltage is expressed in terms of the 
Volterra series of the input voltage; the resulting vector matrix equations 
are successively solved. Since the procedure is similar, the details are 
omitted. 

Two common-collector stages were cascaded using this approach. 
(Fig. 18). The measured values at 120 mA, 10 V with 75 ohm source and 
load impedances were —87 dBm and —112 dBm for the second and 
the third harmonic distortion, respectively. The computed distortion 
values are -88.7 dBm for second and -116.6 dBm for third harmonic 
distortion. Thus, good agreement with experimental result is obtained. 

The cascade formulas are simple, physically meaningful and yield 
rules of thumb to judge the effect of the first stages. The nodal approach 
is more complicated. However, the advantage of the nodal approach 
is that it is general and can be used for an amplifier. For example, a 
cascade of common-emitter and common-collector stages involves five 
nodes; if shunt feedback is used at the input and at the output, the 
same program can be used to analyze this amplifier. (Cascade formulas 
do not take feedback into account.) In general, the nodal approach can 
be extended to study frequency-dependent nonlinear network with 
n nodes, if the nonlinearities are small. 




e i* ©t 

|Y(-v 2 ) 



9(-v 4 ) 



g(v 2 -v 3 ,-v,) g(v 4 -v 5) -v 3 ) 

Fig. 18 — Common-collector — common-collector nonlinear equivalent circuit. 
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VII. ENGINEERING APPLICATIONS 

A few pertinent practical applications of the work are described 
below. These results were either first predicted by the model and then 
verified in the laboratory or first experimentally observed and then 
confirmed by analysis. 

In the initial design of L4 repeater a common-emitter — common- 
emitter — common-collector configuration 28 was used in the power am- 
plifier. The third harmonic modulation performance was not as good as 
desired. This led first to the study of the output common-collector stage. 
As shown in Fig. 19, the increase in source impedance increases the 
distortion of the common-collector stage. Since the preceding common- 
emitter stage output impedance is high, the common-collector per- 
formance was not optimum. Secondly, the preceding common-emitter 
stage was studied because the gain of the common-collector stage is low. 
(see Section VI) As shown in Fig. 20, increase in load impedance beyond 
optimum R L degrades its performance radically. Since the common- 
collector input impedance is high, the common-emitter stage perform- 
ance was not optimum either. Thus, in the redesign work by Ken 
Tantarelli, the common-collector output stage is not being used. 

Another interesting application feature was the improvement in 
modulation performance of the common-emitter stage with increase 
in voltage. As shown in Fig. 8, it is a function of load impedance, and 
at about 150J2, maximum improvement was obtained. 

New coaxial systems are currently being studied to operate at higher 
frequencies. Different configurations have been examined for the output 
stage. The model showed that common-collector and common-base 
performance is poorer with an increase in frequency and thus the use 
of these stages as output stages was questioned (unless transistors with 
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higher /«'s are available). Recently when a new, high-frequency modula- 
tion test set was built, experiments confirmed the prediction. The third 
harmonic coefficient (MS) for a + b — c product was 8 dB poorer at 

36.5 MHz (due to signals at 36.5 MHz, 40.1 MHz, and 43.1 MHz) 
compared to the value at 17.3 MHz (due to signals at 14.5, 15.2, and 

16.6 MHz). The common-emitter configuration modulation performance 
suffered only about one dB degradation. 
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APPENDIX A 



A.l Higher-Dimensional Transforms 

The second-degree case is illustrated as an example. From (4), 

y 2 (t) = [ [ c 2 (t - Tl , t - r a ) II x(r t ) dr { . (40) 

Jo •'0 «=1 

If the system is physically realizable, c 2 (t — t v , t — t 2 ) = 0, for t, > t. 
Hence, the limits can be extended to a> . 

nOO /!» 2 

y 2 (t) = I c 2 (t - Tl , t - t 8 ) II x(T t ) d Ti . (41) 

JO Jo i-1 

Introducing dummy variables t x and t 2 , the two-dimensional trans- 
form is taken 

7 2 (s, ,s 2 )=l I ytfo , t 2 ) exp (-Sjij) exp (-s 2 J 2 ) dt x dU 

Jo •'o 

= / df, / dt 2 I dr t l dr 2 c 2 (t x — T\ , t 2 — t 2 ) 
Jo Jo -Jo Jo 

• II as(r<) dr t exp (—«,«,) exp (-s 2 < 2 ). (42) 

i-1 J 

Substituting <! — tj = wi! , t 2 — t 2 = m 2 , and using the fact that 
c 2 (mi , m 2 ) = for m< < yields 

f%00 /»QO rtCO /»CO 

Y 2 (s, ,s 2 )= I dm, / dm 2 I dr x I dr 2 c 2 {m x , m 2 )x{r ,)x{t 2 ) . 

Jo Jo Jo Jo 

•exp (— s,™,) exp (— 8,t0 exp (— s 2 m 2 ) exp (— s 2 t 2 ) (43) 
= C,(«, , s 2 ) X («,) X (s 2 ). (44) 

A.2 T/ie Output of the Kernels to Si?iusoidal Inputs 

For the second-degree case, consider two sinusoidal signals at fre- 
quencies f a and f b . The input x (r) "equals, 

, x Texp (ju a T ) + exp (-jco n T) "1 f exp (p} b r) + exp (-m,t) 1 
«(r) = _ 2 J + L 2~ J 

(45) 

From (41) 



2/(0 = / dr! / dr 2 c 2 (t — Ti ,t — t 2 ) 

Jq Jq 
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f cxp (p)„T x ) + exp {-pj a T X ) exp ( mtQ + exp ( — ju h r l ) 1 
"L 2 + 2 J 

exp (jo) a T 2 ) + exp (-jco„t 2 ) exp (jo>,,t 2 ) + exp (~m,t 2 ) 1 _ . 

Considering one cross term only, 

,. 30 /too 

/ dr, / dr 2 c 2 (/ — r, , / — r_>) 

•| exp (jo} a T t ) exp 0co 6 T 2 ). (47) 

Substituting ?n, = / — t x , »t 2 = £ — r 2 and carrying out the integration 
yields 

lC 2 (ju a ; p> b ) exp [;(«. + *)<]. (48) 

This term occurs twice as does its complex conjugate. 
Hence, the output due to the a + b term alone is 

y«+»(0 = I C 2 (p} a , ju b ) | cos [(«„ + w 6 )i + a+i]. (49) 

The 2oj a term and its conjugate occur only once in (46); hence, it is 
6 db better. The response of the third harmonic kernel to three sinusoidal 
inputs is similarly treated. 

A.3 Cascade Relations 

For the system shown in Fig. 1, the cascade? formula are given below. 
The cascade relations can be symbolically written as 

Z = F(x) = E{y) = E(D[x\) = (E o D)(x). (50) 

The formula are 

F l (8 1 ) = 8,(0D,(0 (51) 

2 

F,(«, , s 2 ) = E l (s 1 + s 2 )D 2 ( Sl , s 2 ) 4- E 2 ( Sl , s 2 ) II #.(*.) (52) 

F 3 (s l , s 2 , s 3 ) = J5,(sx + s 2 + 8 S )D 8 (*i , s 2 , a 8 ) 

3 

+ 2E 2 (s 1 , s 2 4- s 3 )D 1 (s,)D 2 (s 2 , s 3 ) + # 3 (*i , s 2 , s 3 ) JJ D i( s <)- (53) 

• = i 

A physical interpretation of the formula for cascaded transistors is 
given in Section VI. The procedure for deriving the cascade relation is 
as follows: the output Z(t) of the last stage is expressed in terms of the 
Volterra series of its input. (Only two terms are considered) 



TRANSISTOR DISTORTION ANALYSIS 1021 

Z(t) = f e t (t - r,MTj dr x 

+ \ l *S - n , t - t 2 ) II y(r,) d Ti . (54) 
The output of the first stage y(t) is related to its input by 
V{r) = di(r - o-O-rCo-,) d<r 

^o 

+ P r rf 2 (r - ^ , r - «r 2 ) II sfa) dr, (55) 

•' J t' = l 

Substituting (55) in (54), terms of the same degree are collected; as an 
example, the first second-degree term equals 

f dr e x (t - r) ff d 2 (r - a, ,t- a 2 )f{ .r(er.) da, . (56) 

Taking the two-dimensional transforms yields 

*,(•, +s 2 )D 2 ( Sl ,* 2 ) fl *(*,). (57) 

i=i 

APPENDIX B 

The Nonlinear Parameters 

From (15), 

L = KIbMV cb ). (58) 

A two-dimensional Taylor's series expansion of (58) is taken; i, is 
expressed by K(v 2 ) and v C n = v 3 — v x . Hence, 

i t = g(v 2 , /•:. - ''.) = 0iO' 2 , ih - i>i) 

+ (JoM , v 3 — »,) + foO'2 , t> 3 - <',), (59) 
where 

gM , v 3 - t>j) = a,MoK>2 + M l (v 3 - v t ), (60) 

f7 2 (r 2 , v 3 - v,) = a 2 M K\{v 2 f + m 2 (v 3 - v t f 

+ ai M K 2 (v 2 ) 2 + a&KMfa - Vl ), (61) 
and 

0a(»a , V a - i\) = a 3 Mo(^i) 3 fe) 3 + m 3 {v 3 - v,) 3 
+ a x M 2 K x v 2 (v 3 - Utf + a-SUK\vl{v 3 - Vi) 
+ aAKM* + 2a 2 il/ /C,/v 2 (i- 2 ) 3 + aAtf^) 2 ^ - v,). (62) 
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The avalanche coefficients are given below. 

M = 



\V CB0 ) 



(63) 



*i = (M )' = n % CB) V (Mo) 2 

K " CBOJ 

V CB J-'J-O 



27, 



0*«/ 



Mi 







(64) 




(65) 


3 


f(n - 1) /^.Vl 
L2(F CB ) 1S + W / J 


•(66) 



The coefficients m x , m 2 , and m 3 equal m { = Ic(Mi/Mo)', i = 1, 2, 3, 
where I c is the collector dc bias current. 
The h FE coefficients are given below: 



«i = 



tips i 



K E max + 1 + a log 2 -j— ^— + 2a log e log ■=-* 



<*• - -S7- r^" 2a log J log y^- + log e 

L * C max 



21 c ilpE 

-2a 2 



«3 = 



6 



+ 12 



(« 2 ) ! 



1 



(a,) 3 



i_ 7 C ' («i) (7 C ) h FBmaT 
The collector capacitance coefficients are given by 

7. = KVcb)-* 

7 2 = ^ ikFm)-* 

A; _ z 

73 = 27 (^cs) 3 - 



2a(log e) : 



■]■ 



(67) 

(68) 
(69) 

(70) 

(71) 

(72) 



From (62) for g 3 (v 2 ; v 3 — Vi), <? 3 is obtained by replacing 2?i(s.) for 
v 2 and Ci(Si) — Ai(s { ) for (v 3 — v0; moreover, the kernel must also be 
symmetrical. Since the procedure is the same as for § 2 it is omitted. 
The interaction terms are given below: 



23 = 2a 2 lft Q K\ B x {s,)B 2 (s 2 , s 3 ) 



+ 2ra 2 [Ci(s,) — ^. 1 (s,)][C 2 (s 2 , s 3 ) — A 2 (s 2 , s,)] 
+ «AX, B 2 ( Sl , ±)[CM - AM] 
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+ cLyMJLy B l (s l )[C 2 (s 2 , s 3 ) - A 2 (s 2 , «,)] (73) 



723 = 2 T2 [CM ~ Bi(«i)][C.(«2 , s 3 ) - £ 2 (s 2 , s 3 )] (74) 

where— denotes symmetrical kernel. 

appendix c 

A2436 is an n-p-n silicon transistor with overlay type of construction. 
Its ] T ranges from 800 to 1000 MHz. It is a power transistor with current 
capability of 1 amp and can handle 2.2 watts of power. 

Typical parameter values for transistor type 2436 27 at 120 mA, 
10V are given below: 

I e =0.12 amps. 

r b = 13.6 ohms 

r e = 5200 ohms 

C, = (6)10~ 12 farads 

C 2 = (3.97)10"° farads 

C 3 = (9.2) 10" 12 farads 

Z„ = 50 ohms 

Z L = 50 ohms 

V CB =10 volts 

Vcbo = 350 volts 

n =2 

re, = 0.2165 ohms 

a = 0.38 

h FE mnx = 122 

Ic max = 0.633 amps. 
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